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ABSTRACT 
For  any  system  whose  input-output  relationship  may  be 
described  by  a  ratio  of  polynomials  in  S  (the  complex  number  7*  £Ju?), 
one  of  the  measures  of  the  system  performance  is  its  transient 
response  to  a  step  input.  For  systems  of  order  greater  than  second 
the  d etermination  of  the  parameters  of  this  transient  response 
has  always  been  a  difficult  process.  In  this  thesis  a  new 
procedure  is  developed  which  allows  the  rapid  and  accurate 
determination  of  the  transient  parameters  of  any  system  up  to 
and  including  third  order  with  one  zero.  The  procedure  consists 
of  developing  the  response  equation  in  terms  of  angles  at  the 
singularities  rather  than  the  magnitudes.  The  resulting  equation 
is  utilized  to  develop  expressions  for  times  of  extremals, 
magnitudes  of  extremals,  and  a  form  of  rise  time.  A  simple  but 
accurate  graphical  method  is  demonstrated  for  solving  the  result- 
ing transcendental  equations.  It  is  shown  that  the  expression  for 
times  of  extremals  is  a  function  of  only  two  variables  and  may 
therefore  be  displayed  on  a  two-dimensional  plot.  The  times 
determined  from  this  plot  are  utilized  in  computing  the  magnitudes 
of  extremals.  Various  curves  are  presented  to  provide  a  means  for 
general  analysis. 
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TABLE  OF  SYMBOLS  AND  ABBREVIATIONS 


Definition 


S  The  complex  frequency  GS'F"^^). 

V  The  real  part  of  a  complex  number; 

neper  frequency. 

Co  The  imaginary  part  of  a  complex 

numberj  radiari  frequency. 

tfy  *y^/  The  S-plane  coordinates  of  a  pair  of 

complex  poles. 

opn  The  natural  frequency  (uJ„  ^\7~ct*'  r  ciy*"), 

"g  The  S -plane  coordinate  of  a  real  zero. 

0  The  S -plane  coordinate  of  a  real  pole. 

TT  The  angle  measured  clock  wise  from 

the  negative  real  axis  to  the  vector 
from  the  origin  to  the  complex  pole 

q  The  angle  measured  counter  clock  wise 

from  the  real  axis  to  the  vector  from 
a  real  pole  to  the  complex  pole  TJ+J&Jc/' 

q  The  angle  measured  counter  clock  wise 

from  the  real  axis  to  the  vector  from 
a  real  zero  to  the  complex  pole  ^-/y-W* 

/  Real  time. 

-7-  Nondimensionalized  time  (T'^t). 

j  The  time  to  the  rt"  point  of  zero  slope 

'*  of  the  transient  response. 


The  time  to  the  n"   point  at  which  the 
transient  response  is  equal  to  the 
final  value. 


A(4)  The  time  response, 

Mf/n  Normalized  magnitude  of  the  transient 

response  at  time  ipfi  (tfpLs  X(ip„) //•**>  ho  he). 


1.   Introduction. 

In  recent  years  the  study  of  systems  which  have  an  input- 
output  relationship  (or  a  relationship  between  any  two  variables 
of  a  system)  which  is  described  by  a  differential  equation  has 
increased  tremendously.  Such  relationships  are  found  in  the 
equations  of  motion  of  objects,  in  the  equations  for  chemical 
reactions,  in  thermodynamic  relations,  in  the  equations  for  heat 
and  fluid  flow,  in  electrical  network  equations,  in  automatic 
control  and  servomec nanism  equations,  and  even  in  the  equations 
governing  the  fields  of  industry  and  management »  These  relation- 
ships are  written  in  either  of  two  ways: 

(a)  Given  that  the  variables  of  interest  are  X(t)  and  Y(t) 
which  are  expressed  as  functions  of  some  third  variable 


W(t)  and  its  derivatives s 

(b)  Given  that  the  variables  of  interest  are  X(t)  and  Y(t) 
which  are  expressed  implicitly  in  the  equations 

k,J"%r* --+£x-*- <*&"•* **>' 

In  either  case  the  complex  frequency  (5  =  v^/^?)  may  be  introduced 
and  the  relationship  expressed  ass 

/0    *<±i  -  ^lll^llZll —  *  a° 

The  usual  case  of  interest  is  to  allow  Y(t)  to  be  some  input 
(either  an  initial  condition  or  a  specific  function)  and  to 
investigate  the  response  of  X(t).  A  specific  case  which  is  often 
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mentioned  in  design  specifications  is  the  response  of  X(t)  when 
Y(t)  is  a  unit  step  input  (i.e.,  a  unit  value  input  applied 
instantaneously  at  time  zero),.  Depending  on  the  order  of  the 
relationship  ,  the  value  of M  in  the  equation  (l),  the  determinat- 
ion of  the  response  may  be  straightforward  or  quite  tedious  and 
time  consuming.  For  the  pure  second-order  system  (m=2  and  n=0) 
the  response  has  been  extensively  investigated,  and  general  curves 
and  equations  are  available  in  any  elementary  text  on  servo- 
mechanisms.  For  relationships  with  m>2  and/or  n?0  the  information 
is  not  nearly  so  plentiful.  The  approach  taken  by  most  authors 
is  to  present  the  equation  in  basic  form,  comment  briefly  on  the 
effects  of  the  higher  order  terms,  and  turn  to  methods  of  making 
these  higher  order  terms  negligible  in  order  that  the  response  may 
be  analyzed  in  terms  of  the  more  familiar  second-order  relation- 
ships. Unfortunately,  the  higher  order  terms  cannot  always  be 
made  negligible,  and  the  engineer's  only  recourse  has  been  to  work 
out  each  system  in  detail.  No  method  has  been  available  for  a 
general  analysis  of  the  effects  of  the  higher  order  terms.  In  this 
thesis  a  procedure  is  developed  which  enables  the  analyst  or  designer 
to  solve  in  general  an  equation  of  one  higher  order  in  both 
numerator  and  denominator  (i.e.,  order  three  in  denominator  and 
order  one  in  numerator.) 

Perhaps  the  simplest  way  to  gain  an  appreciation  of  the  problem 
is  to  examine  the  relationship  on  the  S-plane.  Equation  (1)  may 
be  rewritten  as  follows: 


Where  the  Zi  and  Pi  may  have  either  or  both  a  real  part  and  an 
imaginary  part. 

An  example  of  such  a  system  appears  in  Figure  1. 

a 


a-  4-   J^i 


i 


ft* 


D ,_ 
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Figure  lr  An  example  of  S -plane  diagramming 
The  system  depicted  would  be  written  as: 

where  the  S  term  in  the  denominator  represents  Pin,  the  step  input. 

The  usual  nomenclature  assigns  the  term  poles  to  the  P 
values  (represented  by  D  )  and  the  term  zeroes  to  the  Z  values 
(represented  by  o)  with  the  collective  name  of  singularities  assign- 
ed to  all  roots  of  both  the  numerator  and  denominator  polynomials. 
These  singularities  may  occur  anywhere  in  the  S -plane  and  in  any 
combination,  but  for  stable  physical  systems  the  poles  always  occur 
in  the  left  half  plane  (vV£) ,  and  complex  singularities  always 
occur  in  conjugate  pairs. 

The  second-order  approximation  approach  is  to  consider  the 
"dominant"  pair  of  complex  poles  (7^^k*^in  the  example)  as  the 
determining  factor  and  force  all  other  singularities  far  out  in 


the  left  half  plane  so  that  they  may  be  considered  negligible, 

A  generalized  method  of  analysis  for  systems  with  one  pair 
of  complex  poles  and  one  real  singularity  has  been  developed 
by  Patton  and  Abbott  [l]  .     This  method  consists  of  using  a 
scale  factor  to  transform  the  real  singularity  to  the  Ts  -J 
position  on  the  real  axis  and  then  determine  various  response 
parameters  by  entering  a  set  of  generalized  curves  with  the 
location  of  the  scaled  complex  poles  in  relation  to  the  fixed 
singularity.  A  study  of  this  method  was  the  starting  point  in 
the  search  for  a  method  of  generalizing  the  equation  for  a  system 
with  a  pair  of  complex  poles,  one  real  pole,  and  one  real  zero. 


2 „   Attempts  at  Generalizing  the  Equation  for  a  Third  Order 

System  with  One  Zero. 

A  clarification  should  be  made  of  what  is  meant  by  general- 
izing the  equation.  If  the  stated  problem  contains  more  than 
two  variables,  it  is  desirable  to  combine  the  variables  by  means 
of  some  linear  transformation  so  that  a  maximum  of  three  trans- 
formed quantities  remain  as  variables  and  all  other  variables 
are  "buried"  in  the  transformation. 

This  point  is  illustrated  by  means  of  two  examples: 

(1)  The  pure  second-order  response  to  a  unit  step  input 
is  expressed  ass 

v       -  £ r- 

acs;  "  stsi+tK/s  +  v**) 

Assume  that  a  time  (time  of  peak,  rise  time,  settling 
time,  etc.)  is  the  transient  parameter  of  interest. 
There  are  then  three  variables  of  interests  ^/,  <^,f  / 
(k  appears  only  as  a  multiplier).  It  is  possible  to 
display  the  information  desired  on  a  two-dimensional 
graph  by  displaying  %/   on  the  abscissa,  t  on  the 
ordinate,  and  ty  as  a  parameter  for  a  family  of  curves: 
however, the  problem  may  be  simplified  by  noting  that 
the  time  solution  may  be  written  ass 

*w=  s&p-^^P2  *'r^'"  ^*T)] 


Now  let  T-4k/,  and  U-  &/u?J  : 

7T=  fire  /an    (lit) 
This  information  could  be  displayed  on  a  two-dimensional 
graph  with  U  as  abcissa  and  T  as  ordinate.  The  necessity 
for  the  family  of  curves  has  been  precluded  by  "burying" 
the  variable  &)J  in  the  transformation  T=^/. 
(2)  As  a  second  example,  consider  the  response  of  the  pure  third- 
order  system  to  a  unit  step  input: 

j  £_ - 

X  (i)    -  4(J*P)(S*+ZTSS'T&nl) 

There  are  now  four  variables:  TU,  ^/3  /  ,'/**'•  Obviously,  four 
variables  cannot  be  displayed  on  a  two  dimensional  plane. 
Writing  the  response  in  the  time  domain:      A  — , 

Now  let  T=i^/J  t   U=  ^,/and  V=  %-*/:  

*7^)  =  -%  fj   -  --^-A  *"*r-  J**"*'— >'-*"*■-  crtr-dtil 

f*firch>(!<*)   5  &f=  0rcJ**   Cvh) 

This  information  may  be  displayed  on  a  two-dimensional  graph 

by  plotting  U  on  the  abcissa,  T  on  the  ordinate  and  V  as  a 
parameter  for  a  family  curves.  Again,  the  fourth  variable  has 
been  "buried"  in  the  transformation  T=^//. 
It  may  be  noted  that  the  transformation  used  in  the  above  examples 
was  not  the  only  possibility.  In  either  examples,  another 
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possibility  would  have  been  T=Ttj  and  in  the  case  of  the  third 
order  system,  another  possibility  would  have  been  1=Pt.     The 
latter  is  the  transformation  used  by  Patton  and  Abbott  [J ]  who 
showed  that  this  transformation  is  equivalent  to  placing  the  pole 
at  minus  one  on  a  W-plane  defined  by  W*s/p.  These  authors  then 
displayed  the  response  parameters  of  rise  time,  time  of  peak 
overshoot,  magnitude  of  peak  overshoot,  and  settling  time  by 
providing  curves  on  the  W-plane  representing  constant  values  of 
these  parameters.  The  parameter  values  for  a  particular  system 
may  be  determined  by  locating  the  conjugate  roots  on  the  W-plane 
and  reading  the  value  of  the  curve  passing  through  these  points. 
The  response  of  a  third-order  system  with  one  zero  to  a  unit 

step  input  is  expressed  asr 

u        _  _JUJLLZ2— 

ACS)-      £(stP)(<S*'+  ZTcJS  +  t*)^) 

The  number  of  variables  has  now  been  increased  to  fiver  u)g/^%J^ r9  P*f  2. 
This  number  can  be  reduced  to  four  by  means  of  a  transformation 
similar  to  the  type  just  described,  but  a  maximum  of  three  may  be 
displayed  on  a  two-dimensional  curve.  In  order  to  reduce  the 
problem  a  transformation  must  be  found  which  "buries"  two  of  the 
variables.  One  approach  is  made  by  noting  that  the  pure  third 
order  system  was  generalized  by  making  a  transformation  which 
placed  the  real  pole  at  the  minus  one  position  on  a  W-plane. 
It  appears  that  if  a  transformation  could  be  found  which  would 
fix  the  two  real  singularities  at  two  predetermined  points  on  some 


W-plane  where  W  is  a  linear  function  of  S  then  the  only  variables 
would  be  the  angle  and  distance  of  the  conjugate  poles  from  the 
origin  of  this  W  plane.  In  order  for  this  transformation  to  be 
useful  there  must  exist  a  unique  relationship  between  the 
fictitious  time  associated  with  the  W-plane  and  real  time  associated 
with  the  S-plane,  and  likewise  such  a  relationship  must  exist 
between  fictitious  magnitude  and  real  magnitude. 

Although  many  minor  attempts  were  made  to  determine  a  trans- 
formation of  this  type,  only  two  seemed  to  have  any  merit t 
(1)  The  bilinear  transformation. 

It  is  well  known  that  by  means  of  the  bilinear  transfor- 
mation any  three  points  on  the  S-plane  can  be  mapped 
into  any  three  points  on  the  W-plane.  If  it  is  desired 
to  map  the  points  P,  Z,  and  origin  into  the  points  -1, 
-2,  and  the  origin,  this  may  be  achieved  by  the 
transformation? 

which  reduces  tor. 


W  = 


(2)  Transformation  by  scaling  and  translation. 

Another  means  of  placing  the  pole  at  the  minus  one 
point  and  the  zero  at  the  minus  two  point  is  by 
compressing  (or  expanding)  the  scale  until  the  distance 
from  the  pole  to  the  minus  one  point  is  equal  to  the 


distance  from  the  zero  to  the  minus  two  point  and 
then  translating  through  this  distance.  This  may  be 
accomplished  by  the  transformations 

Although  either  of  these  transformations  can  be  shown  to 
be  suitable  for  certain  specific  cases  no  unique  relationships 
could  be  determined  which  would  apply  for  the  general  case. 

No  f  urther  attempts  were  made  at  determining  a  trans- 
formation. Although  no  proof  is  given,  it  will  be  shown  in  the 
next  section  that  it  is  unlikely  that  a  suitable  transformation 
exists . 

One  of  the  major  difficulties  experienced  by  Patton  and  Abbott 
even  in  view  of  the  generalizing  transformation,  was  the  solution 
of  a  transcendental  equation  of  the  form: 

jr,  e         =  S'n  (tai  *  tz) 
In  computing  values  for  the  curves  in  their  thesis  they  were  forced 
to  use  iterative  methods  for  the  solution.  While  studying  their 
work  a  graphical  method  was  discovered  for  solving  this  form 
of  equation,  and  this  discovery  led  to  the  idea  for  a  means  of 
generalizing  the  third  order  system. 


3.   Development  of  the  General  Equation, 

The  response  of  the  third  order  system  with  one  zero  is 


expressed  in  the  S  domain  as: 

AOiJSl 


k  (3*2) 


X(s)=  j(siP)(j>*rZfy6  +  wS)     "     *  Uirt&rJyMfrV-ji+t) 


and  is  expressed  in  the  time  domain  as: 


X<v;- 


toiiMJ. 


P 


/ 


-  ±Al?^L~  $£jw^> f  (v-Wrt *  ^  Co^  ^jj. 

The  time  domain  equation  is  derived  by  simple  residue  methods  and 
algebraic  manipulation.  The  various  values  are  depicted  in 
Figure  2. 


Figure  2:  S -plane  plot  of 


S+  Z 


The  expression  for  X(t)  may  be  simplified  by  rewriting  it  in 
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terms  of  the  angles  ©p,  0z,  and  Tf  ,  The  following  relationships 
are  written  by  inspection  from  Fig,  2: 

21?   -  (Z'Vf  +  uht*' 
S,  n  7T  -     ^/»>n  ;  COS  2T2    ^/t^rn 

S,n<p/oz    wJ/poi  C&J&/9Z    P-Toyffi 

By  trigonometric  identities: 

From  these  relationships:  m<Sjl 

By  further  use  of  the  relationships  above: 


Si"*P- 


O  >  '7 


And :  y 

Therefore  :  -fo/^  +&  Uj<*" 

Eq'n.  I  has  a  number  of  important  properties  which  are  noted  here: 
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(1)  For  stable  systems  the  multiplying  constant  -jJ^fAxs   the 

final  value  determined  bys 

Final  Value  -  Jb*    J  X(s) 
■3-*o 

(2)  If  the  transformation  T»4^  is  made  the  equation  is 
completely  non-dimensionalized|  i.e.,  the  equation  is 
a  function  of  angles  only# 

(3)  The  equation  may  be  normalized  to  unity  by  multiplying 
by  k^-     jrg    j  i.e.  the  final  value  has  been 
normalized  to  unity.  Using  these  properties?      , 

Eq'n  la  is  valid  for  8z#L80  -  T  and  ©p  <   l80#-  f  j  i.e.  for  the 
zero  not  at  the  origin  and  the  pole  in  the  left  half  plane  (it 
is  assumed  that  T£  "%)•     When  ©p  £l8Q -^the  system  is  unstable 
and  the  final  value  theorem  does  not  hold.  When  ©z-l80  -  V       the 
zero  is  at  the  origin  and  the  final  value  is  zero.  These  special 
cases  are  treated  in  detail  in  later  sections. 

Eq'n  la  is  valid  for  any  order  system  up  to  and  including  third 
order  with  one  zero.  Any  combination  may  be  achieved  by  allow- 
ing either  or  both  of  the  angles  0p  and  6z  to  go  to  zero  which  is 
equivalent  to  locating  the  singularity  at  infinity. 
This  equation  may  now  be  manipulated  to  determine  various  parameters 
of  the  time  response. 
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U.   Development  of  the  Equations  for  the  Parameters  of  the  Time 

Response. 

For  a  given  system  the  slope  of  the  response  curve  at  time 
zero  can  be  quickly  computed  along  with  the  final  value  of  the 
response.  Since  it  is  also  known  that  any  oscillatory  motion 
is  of  the  damped  sinusoid  form  it  is  necessary  to  know  only  three 
other  parameters  to  be  able  to  sketch  quite  accurately  the 
complete  transient  response.  These  parameters  ares 

(1)  tc-time  at  which  the  response  crosses  the  final  value. 

(2)  tp-time  of  extremals  (points  where  slope  of  time  response 
is  zero.) 

(3)  Mpt-normalized  magnitude  of  the  response  at  the  extremals. 

(a)  Development  of  equation  for  the  time  of  extremals  (tp): 
In  order  to  determine  this  value  the  expression  for 
X(t)  is  differentiated  and  equated  to  zero? 


equating  to  zero: 


sos  -i#t't/h**p  . 

13 


There  is  no  "closed"  method  for  solving  this  equation  for  tp> 
but  a  simple  graphical  method  exists  and  will  be  explained 
presently.  Note  that  the  times  of  the  extremals  are  independent 
of  the  angle  If, 

(b)  Developments  of  the  equation  for  normalized  magnitude  of 
extremals  (iMpt)r  To  determine  this  value  the  equation  for 
X(t)  is  entered  with  the  results  of  part  (a): 

fa  -    yWii2   =  $1     +    ±n_%>  St„C*j  -*f)       -  (P-JW     -y> 

"r  r/twsr5      Z^    >s,*t  s*nC&i.tr)  s 

_  s,n(&fi?)_  emrr^.„  (*t/Jfi6>i 'Spiral 
s^rSi^ceztir)  ^j 

Substituting  from  Eq'n  II  for  the  quantity  in  braces s 


Values  for  solution  of  this  equation  are  determined  graphically 
at  the  same  time  Eq'n  II  is  solved, 

(c)  Development  of  the  equation  for  the  times  when  the  time 
response  crosses  the  final  value  (tc)t 
This  time  occurs  when  the  last  two  terms  of  Eq'n  I  are 
equated  to  zero: 

S-nTT  S>r>(&ztr)  S>iTT  S>n(t9iiS) 


S»oT  S-r>C&ztr)  S>njr  S,n(<9iia) 

0  <;>n(6>ptir} 


Again  tc  may  be  determined  graphically  on  the  same  plot  used  to 
solve  Eq'ns  II  and  III. 

The  notation  of  the  equations  is  simplified  by  making  the 
transformations  as  before s 


/*-  ""^i 


Hi 


The  transformed  equations  then  becomes 


Uf>.  n* :    S'»irS"><£L£!L  e -  *%** m 


-  St*  (Tf+<9i~fy>  tr>J< 


S,s>  (&fiT) 


*    £>*  (Tc-t&£  -<?/>  i  T) 
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5.  Graphical  Solution  of  the  Equations. 

Equations  Ila  and  IVa  may  be  solved  graphically  to  an  accuracy 
of  approximately  one  percent  by  following  the  step  by  step 
procedure  outlined  below.  Equation  Ilia  may  be  solved  directly  after 
the  solution  of  Equation  Ila  for  Tp. 

(a)  Graphical  solution  is  performed  on  semi -log  paper  with 
values  of  T  plotted  on  the  linear  axis  and  values  of  the 
functions  of  T  plotted  on  the  logarithmic  axis. 

(b)  Compute  or  measure  the  angles  0p,  9z,  andT  on  the  S -plane 
diagram  of  the  problem  to  be  solved. 

(c)  Compute  the  following  values: 

Sin  (0z=6p) 

Sin  0p/sin  (0p+7O 

(d)  If  Sin(0z-0p) >  0   Construct  the  function  Sin  T. 
If  Sin(6z-ep)<  O    Construct  the  function  Sin(-T)* 
(Note:  This  function  is  easily  constructed  by  plotting 

the  values  of  Sin  T  for  T=  10°,  20°,  30°,  U5°, 
60°,  90°  and  appropriate  multiples.  For  T  <  10° 
the  function  is  approximately  a  straight  line.) 

(e)  Construct  the  function  Sin(6z-9p)  &  • 

(Notes  On  the  semi-log  paper  this  function  is  a  straight 

line  with  intercept  Sin(0z-0p)  and  slope  of  -  JJ^$2   .) 

(f )  Locate  the  point  on  the  straight  line  where  the  directed 

horizontal  distance  to  the  sine  function  is  equal  to 
(0z-0p).  The  value  of  T  at  this  point  is  Tp. 
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(g)  Construct  the  function  -* ^in  (ep+r)  ° 

(Note?  If  S'nfffferfriis   of  the  opposite  sign  to 

Sin(6z-6p)  it  will  be  necessary  to  make  a  second 
plot  on  semi -log  paper  following  the  rules  laid 
down  in  (d)  above.) 
(h)  Locate  the  point  on  the  second  straight  line  where  the 
directed  horizontal  distance  to  the  sine  function  is 
equal  to  (0z-0p  +  r).  The  value  of  T  at  this  point  is  Tc. 
(i)  Compute  the  value  of  Mpt  by  use  of  Equation  Ila  and  the 

value  of  Tp. 
(j)  Use  the  inverse  transformations s 


x^>-  ^r  Mp*  =  *&*/>+ 


to  obtain  true  parameter  values, 
(k)  Repeat  steps  (f),  (h),  and  (i)  for  as  many  values  as 

necessary  to  obtain  a  complete  description  of  the  transient 
response. 
The  crucial  step  in  the  graphical  solution  just  outlined  is  the 
determination  of  the  sign  placed  on  the  values  of  the  logarithmic 
axis.  Since  the  construction  is  performed ^on  semi -log  paper  it  is 
impossible  to  show  both  positive  and  negative  values  on  the  same 
plotj  consequently,  only  alternate  half -cycles  of  SinT  may  be 
plotted.  If  9z  <  0p  the  value  of  sin(6z-6p)  is  negative  and  the 
logarithmic  axis  must  be  labelled  with  negative  values;  therefore, 
the  negative  portion  of  Sint  must  be  plotted.  Likewise,  if  0z  >  6p 
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the  value  of  Sin(9z-9p)  is  positive  and  positive  portions 

of  SinT  must  be  plotted.  Similar  rules  apply  to  the  solution 

of  Equation  IVa. 

The  procedure  is  illustrated  by  means  of  the  following 

example  t 


Given  a  system  for  which: 


e+  =  *0*  cjjrl 


The  computation  proceeds  as  follows  (see  Fig.  3): 

(1)  Construct  the  function  SinT  on  semi -log  paper. 
(Note:  9z-9p  <  O  j  therefore,  plot  the  negative 

half -cycles ) . 

(2)  Construct  the  line  Sin(9z-9p)  e 

(Notes  9z-6p<0;  therefore,  Sin(9z-9p)<  o) 

(3)  The  equation  to  be  solved  is: 

Locate  the  points  where  the  horizontal  distance 

between  the  straight  line  and  the  sine  function  plot 

is  U0°  (straight  line  point  to  the  right  of  the  sine 

function  point) .  The  T  axis  value  of  the  points  on 

the  straight  line  are  the  values  Tp: 
Tp,  -  4.i    radio** 
7>j  i    (,$2    rod,*** 

~Tp3  -     I C,  3     reo/>e>ri£ 

Note  that  these  values  do  not  depend  on  the  value  of 
T  but  only  on  the  relative  location  of  the  complex 
poles  and  the  real  singularities. 
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Fig. 3:  Cfcmiputational  Diagram  for  the 
System:  Qz*hO;  ep-So";  4"«44? 
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(U)  Construct  the  line;  . 

s,n  r  s>r>cP2tr) 

Read  off  the  values  of  this  line  at  the  times  Tpns 
(Call  the  values  Cn.) 

C-l  -  .213 

C2  -  .165 

C3  =  .113 

(Note:  These  values  are  positive  since  the  logarithmic 

axis  was  assumed  positive  for  construction  of 

the  line.) 

(5)  Read  the  values  from  the  sine  function  at  the 
correspondence  points  found  in  step  (3)s  (Call 
these  values  Bn.) 

B^--  .31 
B2—  .19 
B3—  .103 
These  are  the  values  of  Sin(Tpn+8z~©p). 

(6)  Add  the  distance  (0z-8p+9*)"UU°  to  each  value  of  Tpn 
and  read  the  corresponding  sine  values  (call  the 
value  An#) 

A1—  1.0 
Ag—  1.0 

A3=    Sinii33°=sin730ss.9^ 
These  are  the  values  Sin(Tpn+©z-©p+V). 
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(7)  Compute  the   normalized  magnitudes  of  the  extremals: 


=  1    +    Cn  [3.Sg  3n-/)nJ 


Mpt/»  .9766 
Mpt2=  .730 
Mpt£  =1.071$  . 

(8)  Construct  the  line  , 

-  T/S.U 


<s,n  (ep  +  r) 
=•   -Z.3  & 

(9)  The  equation  to  be  solved  is: 


-Z.ze'^"  =  J<„  (Tf44-°) 


Proceed  as  in  step  (3)  to  find  these  points  of  correspondence. 
The  value  of  T  at  these  points  corresponds  to  the  times  when 

the  response  crosses  the  final  value  (Tc): 

Tc,  «  9.1  radians 
Tc*.«ll„6  radians 

This  completes  the  values  necessary. 


Figure  U:  Step  response  for  &z-  4-0° j  fy-  fr\  T-7+";  &J*i 
In  order  to  demonstrate  the  accuracy  of  the  method,  this 
system  was  simulated  on  the  CDCI6OI4.  digital  computer  using  a 


21 


method  of  solution  completely  independent  of  the  results  obtained 
here.  The  computer  response  is  shown  in  Figure  5  (Notet^d  on  the 
computer  simulation  was  2.0).  For  purposes  of  comparison,  the 
computed  values  and  the  computer  response  values  are  tabulated: 


Parameter 

Computed 

Computer  Response 

Difference 

Tp, 

lul 

U.3 

6.8g 

Tp* 

6.82 

6.9 

2.S* 

Tp,, 

10.3 

10.3 

.96* 

Mpt, 

.9766 

.968 

1.2* 

Mpt* 

.730 

.728 

.55* 

Mpta 

1.0715 

1.076 

.51* 

Tc, 

9.1 

9.1 

1.1* 

Tc* 

11.6 

11.5 

1.75* 

(Noter  Computer  response  values  were  interpolated  from 
the  computer  readout  tabulation  given  in  0.1 
second  steps) 
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Fig. 5";  Computer  Response  for  the 
System:  9z-40;Gp-80;T-44, 
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6.   Adaptation  of  the  Equations  to  Other  Singularity  Combinations. 

It  has  been  stated  that  the  equations  developed  are  general 
in  that  they  are  correct  for  pure  second-order  systems,  second- 
order  systems  with  one  zero,  pure  third  order  systems,  and  third 
order  systems  with  one  zero.  This  property  is  developed  in 
more  detail  in  this  and  later  sections.  Also,  although  an 
important  value  of  the  work  just  described  is  the  method  of  solv- 
ing the  transcendental  equations  to  obtain  various  parameters  of 
the  time  response,  an  equally  important  value  is  that  the  means 
have  been  provided  to  develop  curves  and  relationships  which 
allow  the  analyst  and/or  designer  to  observe  the  effects  of  the 
real  singularities  at  a  glance*  This  point  is  also  developed  in 
detail  in  later  sections. 

To  make  equations  I  through  IV  applicable  to  the  particular 

system  of  concern,  it  is  necessary  to  modify  the  multiplying 

constant   £  ^^/i^nf\    •  It  is  more  correct  to  define  the 

multiplying  constant  as  the  "Final  Value"  of  the  system.  The 

final  value  theorem  for  stable  systems  states: 
M*r>      id)     -  J!Ctr>  3  Pes) 

The  final  values  of  the  various  systems  are  tabulated  in  Table 

1  along  with  the  values  of  the  angles  ©z  and  ©p. 


Svstem  Order' 

Zero 

Final  Value 

©z 

©D 

Two 

None 

t/Vn7- 

Zero 

Zero 

Two 

One 

£2/*>/.** 

Variable 

Zero 

Three 

None 

£  /*>»*? 

Zero 

Variable 

Three 

One 

#B /*>»'*? 

Variable 

Variable 

Table  1:  Modifying  Relationships  for  Equations  I  through  TV. 
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As  an  example  equations  I  through  IV  (modified)  are  listed 
for  order  two  with  no  zeros  (pure  second  order): 

£$n.  X  (mod) ."  Xn )  -  jg>  l  J  'jfcr  C  J/"  (*>Ji  *  r)j 


o 


'/>.*<*»**'):    °   -  J'"  &*/>   =^  Mb»z  n'rr 


-  Mr) "IT  -. 

?f* & (*»</)''  *r^  *'e  Cc* n7T 

These  equations  are  well  known  for  the  second  order  system 
although  the  form  may  be  different. 
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7.   Analysis  of  the  second  order  system  with  one  zero. 
Following  Table  1  of  the  last  section:     / 

£f».*tmJ>:   tyh  [1-   -£*^  j,*^  f^rtf 

which  are  the  governing  equations  of  the  second  order  system  with 

one  zero.  It  may  be  noted  that  each  of  the  time  response  parameters 

may  be  expressed  as  the  pure  second  order  parameter  with  a 

correction  factor.  In  the  case  of  the  times  of  extremals  and 

times  when  response  crosses  the  final  value,  the  correction  is  a 

linear  one;  i.e.,  the  time  is  decreased  by  an  amount  **%^/.  The 

correction  factor  for  the  magnitude  of  extremals  is  a  multiplicative 

nonlinear  factor.  The  magnitude  of  extremals  may  be  written: 

X(tp)  =  Final  Value  (1  +  %   overshoot)       . 

-(celr)*1*' 
andt  %   overshoot  (second  order)  =  -  COS (/)7f)€ 

andt  %   overshoot  (second  order  with  one  zero)  *  -  COS  fflf  ( "   £ 

Vv         /I,         j    /  1   *>«r  '(«*)** 

the  correction  factor  «["j'.ir   e**ft'*r"J 

Figure  6  is  a  plot  of  Mpt  of  a  second  order  system  with  one 

zero  versus  COS  7f  for  various  values  of  Qz*     The  curve  for  9z=0  is 

a  plot  of  the  percent  overshoot  of  a  pure  second  order  system. 

The  system  fails  when  ©z=77"-T  j  i.e.  when  the  zero  is  at  the 

origin  of  the  S -plane.  At  this  point  the  final  value  is  zero  but 
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computed  percent  overshoot  becomes  infinite  resulting  in  an 
undefined  value.  It  is  necessary  to  return  to  magnitudes  to 
interpret  this  special  case: 

le.  My    G>2-  TT-  T  ;    lb-  u>*    ',    "B-O 

for  the  curves  of  Figure  6^n=l  . 

.   y   j   _  j£_     -To-lr 

The  equations  obtained  here  are  quite  simple  to  use  and  no 
further  analysis  of  the  system  will  be  undertaken  in  this  thesis. 
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Fig. 6:  Magnitude  of  First  Extremal  of  a 
Second  Order  System  with  One  Zero 
(Final  value  normalized  to  unity). 
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8.  Analysis  of  the  pure  third  order  system. 

Following  Table  1  of  Section  6:        .     ,    /  . 

ef» .  *  w; :  x„,  =  jfc  {i  -  ftg-Jc     r   -  *g#&fZ»/^ 

"  \        -rr  "  ^'//^^y»  J  \ 

tffn.ll   (we/)'.      -S.^fye  r    =     S>*  (uJc/Jp-Pf) 

It  is  obvious  that  no  simple  correction  factor  exists  for  the 
parameters  as  it  did  in  the  case  of  the  second  order  system  with 
one  zero.  The  methods  developed  in  section  5  have  been  utilized 
to  develop  curves  for  the  parameters  of  the  time  response,  Figure  7 
is  a  plot  of  Tp  versus  0p  (Note  that  Tp  is  independent  of  the 
angle  If  ) .  The  curve  has  been  extended  to  provide  times  for  the 
occurence  of  the  first  three  extremals  (points  of  zero  slope). 
It  is  interesting  to  note  that  the  time  correction  is  practically 
linear  to  0p=7O°  and  may  be  approximately  expressed  as  Tp=n77"V'  6p  up 
to  that  point.  It  will  be  shown  presently  that  no  overshoot  can 
exist  for  9p^90°,  For  purposes  of  comparison  the  times  of  extremals 
are  shown  for  a  second  order  system  with  one  zero  (dashed  lines). 
Figure  8(a)  is  a  plot  of  the  normalized  magnitude  of  the  first 
extremal  of  the  response  as  a  function  of  the  cosine  of  7T    for 
various  values  of  0p  (final  value=l).  Figures  8(b)  and  8(c)  are 
the  same  type  of  plot,  but  for  the  normalized  magnitudes  of  the 
second  and  third  extremals.  The  curve  entitled  6p=0  is  the 
familiar  second  order  curve.  An  important  point  should  be  noted  here. 
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For  a  second  order  system  the  first  overshoot  is  always  the 
greatest;  but  this  is  not  necessarily  the  case  for  a  third  order 
system.  For  0p>76°,  and  dependent  on  the  value  of  7  ,  the  first 
"overshoot"  may  actually  be  less  than  the  final  value  (in  such 
a  case  the  term  "overshoot"  is  misleading);  for  ep>U5°,  and 
dependent  on  the  value  of  Tf,  the  second  "overshoot"  may  be 
greater  than  the  first.  These  effects  are  readily  apparant  upon 
examination  of  Figures  8(a)  and  8(c).  An  interesting  point  is  the 
case  of  @p>90°$  on  the  semi-log  paper  used  for  computational 
purposes  the  line  defined  by  the  following  equation: 

is  a  line  of  positive  slope  and  can  never  intersect  the  curve 
defined  by  the  equation: 

therefore,  in  accordance  with  the  methods  developed  in  Section  $, 
no  time  of  overshoot  can  occur.  This  implies  that  for  6p>90°, 
regardless  of  the  value  of  If  ,  the  time  response  is  monotonically 
increasing  and  never  exceeds  the  final  value.  Figure  9  is  a  plot 
of  tc  for  the  third  order  system  (the  rise  time,  tc,  is  defined 
to  be  the  first  time  when  the  response  is  exactly  equal  to  the 
final  value).  For  values  of  ©p>76°  it  is  possible,  for  certain 
values  of  V ,   that  the  response  does  not  cross  the  final  value  until 
after  the  first  peak.  It  should  be  remarked  here  that  the  curves 
presented  in  Figures  7-9  are  intended  to  serve  as  a  means  of  general 
analysis  and  not  necessarily  as  computational  aids.  The  analyst 
or  designer  interested  in  exact  values  should  perform  the  calculations 
as  outlined  in  Section  5. 
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An  interesting  point  in  the  study  of  the  third  order 
system  is  the  value  0p=76°.  This  value  appears  to  be  the  optimum 
setting  to  minimize  fluctuations  around  the  final  value  regardless 
of  the  value  of  If  %   the  response  always  reaches  the  final  value  on 
first  rise  and  the  first  undershoot  is  minimum  for  all  JT  • 
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Fig. 9:  Rise  Tijne  of  Third 
Order  Systems, 
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9.   Analysis  of  the  Third  Order  System  with  One  Zero. 

The  character  of  the  response  of  the  third  order  system 
with  one  zero  is  best  analyzed  by  plotting  the  values  of  Tp  as 
shown  in  section  5,   this  information  may  be  depicted  on  a  two- 
dimensional  plot  since  Tp  is  dependent  only  on  the  values  of  0z 
and  9p.  This  plot  may  be  constructed  rather  quickly  and  easily  with 
Equation  II  by  fixing  0p  (which  fixes  the  slope  of  the  straight 
line  plot  on  the  semi-log  paper)  and  determining  Tp  for  various 
values  of  9z.  Curves  of  this  form  are  displayed  in  Figures  10(a), 
10(b)  and  10(c)  for  Tp  corresponding  to  the  first  three  extremals 
of  the  response.  If  the  values  of  0p  and  6z  are  known  these 
curves  may  be  entered  and  the  value  of  Tp  for  each  of  the  first 
three  extremals  may  be  determined.  Once  this  value  is  known,  the 
magnitudes  of  the  peaks  may  be  quickly  computed  from  Equation  III. 
Again,  it  is  emphasized  that  the  accuracy  on  the  curves  is  only 
about  *  5%»     If  accuracy  is  required,  the  computational  methods 
of  section  £  should  be  utilized, 

A  great  deal  of  knowledge  of  the  behavior  of  the  system 
response  can  be  determined  from  a  study  of  the  curves  of  Figure  10. 
This  information  is  summarized  in  part  in  Figure  11.  The  coordinate 
plane  used  in  Figure  10  is  depicted  again  in  Figure  11,  but  the 
Tp  curves  are  deleted  and  the  plane  is  divided  into  regions  as 
follows  s 
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Region  I    <    ' 

t    &z  >  */° 

Region  II   /  *f<<ie0 

Region  III  f  <V  >*** 

Region  IV   f  */>  >?"" 

These  regions  are  formed  by  the  line  6p=90°,  which  is  the  angle 
at  which  the  real  pole  moves  "inside"  the  complex  poles,  and  the 
line  0p=0z  which  is  the  "second  order"  line.  The  important  factors 
of  these  divided  zones  are  as  follows: 

(1)  6p<90°  implies  a  negative  slope  to  the  line  which  is  used 
for  computation  as  described  in  Section  5.  9p*90°  implies 
a  positive  slope. 

(2)  0p?9z  implies  that  Sin(0z-0p)<0,  and  Sin(0z-0p)  is  the 
y-axis  intercept  of  the  line  described  above.  Sp^Sz  implies 
that  Sin(0z~0p)  >0»     If  Sin(0z-0p)<C?then  the  negative  half 
cycles  are  plotted  on  the  computational  diagram,  i.e.  the 
values  of  the  Sin  T  for  T»  TT-k  Z770  3TT  J*  4*T, jW/^etc  . 

If  Sin(0z-0p)^0  then  the  positive  half  cycles  are  plotted 
i.e.  the  values  of  Sin  x  for  T  ■  o  Jo7r32r/37%eic.   (See 
Section  5>)» 

(3)  If  the  slope  of  the  computational  line  Sin(9z-©p)£>  ' 

is  negative  then  the  line  must  intersect  all  the  half  cycles 
of  Sin  T  and  an  infinite  number  of  extremal  points  exist.  If 
the  slope  is  positive  then  it  intersects  the  function  Sin  T 
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at  only  a  finite  number  of  points  and  a  finite  number  of 
extremals  exist  (this  finite  number  may  be  zero). 

(U)  ©p->02  implies  that  the  value  T+©z-0p  lies  to  the  left  of 
the  value  T  on  the  computational  diagram.  If  this 
condition  is  coupled  with  the  condition  that  the  slope  of 
the  computational  line  is  positive  then  the  only  manner 
in  which  the  equality  expressed  by  Equation  II  can  exist 
is  for  the  computational  line  to  "pierce"  the  curve  of 
Sin  T  (the  Sin  T  curve  is  concave  downward;  therefore, 
"piercing"  means  to  intersect  on  the  positive  slope  side  and 
go  "inside"  the  curve.)  Any  line  which  "pierces"  must 
"emerge";  therefore,  the  extremals  always  occur  in  pairs. 
If  the  first  extremal  is  called  an  "overshoot",  then  it 
can  be  stated  that  every  "overshoot"  is  accompanied  by 
an  "undershoot";  however,  it  is  possible  that  the  time  of 
the  "overshoot"  and  time  of  the  "undershoot"  may  coincide. 

(5)  0p<8z  implies  that  the  value  T+0z-0p  lies  to  the  right 
of  the  value  T  on  the  computational  diagram.  If  this 
condition  is  coupled  with  the  condition  that  the  slope 
of  the  computational  line  is  positive  then  the  equality 
expressed  by  Equation  III  can  exist  for  a  "piercing" 
condition,  or  it  can  exist  for  the  condition  in  which  the 
computational  line  approaches  just  to  the  horizontal 
distance  Op-Oz  from  the  curve  Sin  T  and  no  closer.  If  the 
latter  condition  occurs  then  an  "overshoot"  and  an  "undershoot" 
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occur  at  the  same  time,  defined  by  the  particular  value 
of  T,  and  no  extremals  will  occur  after  that  time; the 
slope  of  the  time  response  goes  to  zero,  but  does  not 
change  sign.  The  first  extremal  always  exists,  but 
extremals  after  that  time  occur  in  pairs. 
These  factors  are  of  the  utmost  importance,  and  if  they  are 
not  thoroughly  understood  the  reader  should  review  the  procedures 
outlined  in  Section  5  for  computing  Tp»  These  factors  are  graphically 
illustrated  in  Figures  12,  and  13 •  Figure  12  is  a  sample  computational 
diagram  for  0p<©z.  The  values  of  Sin  Tfor  positive  half  cycles  are 
drawn.  One  sample  computational  line  has  been  constructed  for  the 
case  ©p<90°  (0p=73°  for  the  example).  The  points  where  Equation  II 
is  satisfied  are  depicted  by  short  vertical  lines  on  the  computational 
line.  Four  sample  computational  lines  have  been  drawn  for  the 
case  0p*9O°  (6p=108°  for  the  example).  Each  line  has  been  drawn  for 
the  same  0p  but  a  different  value  of  0z-0p.  For  line  "A"  the 
extremals  exist  as  overshoot-undershoot-overshoot  and  no  further 
extremals.  For  line  "B"  the  first  overshoot  occurs,  but  the  times 
of  the  first  undershoot  and  the  second  overshoot  coincide,  and  no 
further  extremals  exist.  For  line  "C"  only  the  first  overshoot  exists 
with  no  further  extremals.  For  line  "D"  the  value  of  6z-0p  has  reached 
its  maximum  of  72°  (0z  maximum  is  180°  and  0p  is  108°),  and  then 
the  first  overshoot  occurs  at  Tp=0j  no  further  extremals  exist. 
Figure  13  is  a  sample  computational  diagram  for  9p*0z.  The  values 
of  Sin T for  negative  half  cycles  are  drawn.  One  sample  computational 


line  has  been  constructed  for  the  case  ©p<90  (9p=73°  for  the 
example).  Four  sample  computational  lines  have  been  drawn  for 
the  case  6p>90°  (0p=lO8°  for  the  example).  For  line  "A"  two 
overshoots  and  two  undershoots  exist.  For  line  "B"  one  overshoot 
and  undershoot  exist.  For  line  "C"  the  times  of  overshoot  and 
undershoot  coincide.  For  line  "D"  no  extremals  exist.  The 
properties  of  the  regions  on  Figure  13  may  now  be  stated  in  terms 
of  these  factors t 

(A)  Region  I,  Referencing  (1),  (2),  (3),  Line  E-Figure 
12:  All  extremals  exist,  slope  of  the  computational 
line  is  negative,  and  positive  half  cycles  are 
plotted  on  the  computational  diagram. 

(B)  Region  II,  Referencing  (l),  (2),  (3),  Line  E-Figure 
13?  All  extremals  exist,  slope  of  the  computational 
line  is  negative,  and  negative  half  cycles  are 
plotted  on  the  computational  diagram. 

(C)  Region  III,  Referencing  (1),  (2),  (3),  (U),  Lines 
A,B,C,D  -  Figure  13:  A  finite  even  number  of 
extremals  exist,  slope  of  the  computational  line 
is  positive,  negative  half  cycles  are  plotted  on 
the  computational  diagram.  In  region  IIIA,  the 
computational  line  will  always  appear  as  line  D, 
i.e.,  there  can  be  no  overshoot  in  this  region. 
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Fig. 12:  Sample  Computational 
Diagram(©p<£z) • 
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Fig.13:  Sample  Computational 
Diagram(0p?9z) » 
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(D)  Region  IV,  Referencing  (1),  (2),  (3),  (U),  Lines 
A,B,C,D  -  Figure  12:  A  finite  odd  number  of 
extremals  exist,  slope  of  the  computational  line 
is  positive,  positive  half  cycles  are  plotted 
on  the  computational  diagram.  In  region  IVA, 
the  computational  line  will  always  appear  as 
line  C  or  D,  i.e.,  there  can  be  only  one  over- 
shoot and  no  further  extremals  in  this  region. 
Regions  IIIB  and  IVB  are  of  special  interest  because  they 
bound  the  line  0z=0p.  Along  this  line  the  pole  and  zero  of  the 
system  cancel  each  other  and  the  system  acts  as  a  pure  second 
order  system.  An  unsuccessful  attempt  at  cancellation  would 
cause  the  system  to  operate  in  regions  IIIB  or  IVB.  Consider 
first  the  region  IIIB  and  refer  to  Figure  13  and  the  computational 
lines  for  0p^9O°.  Assume  that  9p=108  and  0z  is  to  be  varied  from 
0z=lO8°  to  &z=0°  which  will  be  a  movement  along  the  horizontal  line 
0p=lO8°  through  regions  IIIB  and  IIIA.  Obviously  when  6z=0p, 
Sin  (0z-0p)=O  and  all  extremals  exist  and  occur  at  multiples  of  7T 
(this  is  the  second  order  system).  As  0z  decreases  the  computational 
line  translates  upward  rapidly;  when  0z=lO6.3°  |_Sin(0z-9p)=  -•03J 
the  point  has  been  reached  where  the  second  overshoot  and  second 
undershoot  occur  at  the  same  time  (Tp=10.6rad).  When  0z=96.5° 
(Sin  (Oz-Op)=  -»2J  the  point  has  been  reached  where  the  first 
overshoot  and  first  undershoot  occur  at  the  same  time  (Tp=li.5rad)f 
this  point  is  the  dividing  line  between  region  IIIB  and  region  IIIA. 
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For  ©z<°6.5°  no  overshoot  can  exist.  Consider  next  the  region 
IVB  and  refer  to  Figure  12  and  the  computational  lines  for 
6p>90°.  Again  assume  0p=lO8  and  0z  is  to  be  varied  from 
6z=108°  to  0z=l8O°  which  will  be  a  movement  along  the  horizontal 
line  0p=lO8°  through  regions  IVB  and  IVA.  As  0z  increases  the 
computational  line  translates  upward;  when  ©z=113.7  [Sin  (0z-0p)= 
0.1J  the  point  has  been  reached  where  the  first  undershoot  and 
second  overshoot  occur  at  the  same  time(Tp=7.35  rad ) j  this  point 
is  the  dividing  line  between  regions  IVB  and  IVA.  For  0z>113.7° 
only  one  overshoot  can  occur.  When  0z=l8O  [Sin  (0z-0p)=  .95  J 
The  overshoot  occurs  at  its  minimum  time  of  Tp=1.57»  It  is 
possible  to  continue  the  reasoning  in  this  fashion  for  various 
values  of  0p  and  thus  summarize  the  properties  of  regions  IIIB 
and  IVB.  This  summary  is  presented  in  Table  II. 

The  next  factor  of  interest  for  the  third  order  system  with 
one  zero  is  the  normalized  magnitude  of  the  extremals.  This 
information  cannot  be  plotted  on  a  two-dimensional  plot  since 
Mpt  is  a  function  of  0p,  0z,  and  T  •  It  is  not  difficult  however 
to  predict  the  order  of  magnitude  of  the  overshoot.  Along  the 
line  0p=0z  the  system  is  second  order  and  the  normalized  magnitude 
of  the  overshoot  can  be  read  from  the  second  order  curve.  Along 
the  axis  0p=O°  the  system  operates  as  second  order  with  one  zero 
and  the  normalized  magnitude  of  the  overshoot  can  be  read  from 
Figure  6.  Along  the  axis  0z=O°  the  system  operates  as  pure  third 
order  and  the  normalized  magnitude  of  overshoot  can  be  read  from 
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Figure  8(a),  It  is  also  possible  to  plot  the  line  along 
which  Mpt  is  unity.  This  is  accomplished  by  returning  to 
Eq'n  Ilia  and  equating  to  final  values 

Cos  (Tf1<?Z)    -    C*<>  (Tpi0&  +2r) 
W/><c£    j  j    sa /  j  /re. */  W  ««  n  : 

This  curve  can  be  plotted  by  selecting  a  particular  "Jf  and  plotting 

points  on  the  ~Tf>  curves  by  the  equation; 

9g*  2nV  -TT-  Tp 
Curves  have  been  plotted  on  Figures  10a  and  10c,  for  selected 

values  of  T*>  depicting  the  curve  along  which  Mpt=l  (first 

overshoot  in  Figure  10a  and  second  overshoot  in  Figure  10c). 

Note  that  these  curves  end  at  0z=©p=l8O°  -  T*  which  represents 

the  origin  of  the  S-plane.  Beyond  this  point  the  pole  and  zero 

of  the  system  fall  in  the  right  half  plane,  and  the  system  is 

unstable.  A  limiting  situation  in  determining  magnitude  of  the 

extremals  is  encountered  by  letting  0z=l8O-7*(zero  at  the  origin). 

When  this  occurs  Sin  (©z  +7*)=Sin  180°  =0  and  the  expression  for 

Mpt  becomes  infinite,  but  the  final  value  y£^fi\   goes  to  zero  so 

that  the  expression  is  meaningless.  To  evaluate  the  overshoot 

the  value  Z  is  brought  inside  the  parentheses  in  Equation  III? 


h9 


"  (if)    <*„V  £  ^nTS.rthriiT)  \J>,*(ffir)  i  r        *     /       -u 

/Zr    J  he.     Zero    a  J    /*C     Or^'^  \ 


This  equation  is  good  along  the  line  9z=l80°-7'until  the  point 
0p=l8O°-7"  is  reached.     Bringing  P  inside  the  parentheses.* 


«Vi 


j^jM**^ 


*>„*  e 


0ut     -for  Ggs.  1%0-V  y     E9n'  &   £><£con>e*: 


Si  A  (9(0  4  T 


^^r 


Wr% 


which  is  the  expression  obtained  for  a  pure  second  order  system. 
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Another  limiting  expression  occurs  when  9p=l80-  JT 
(pole  at  the  origin).  When  this  occurs  Sin  (©p+P*  )=Sinl80o»0, 
but  the  final  value  y^rzj   becomes  infinite  so  that  the  expression 
is  meaningless.  To  evaluate  Mpt  the  value  P  is  brought  inside 
the  parentheses  as  before:       /  > 

•St*k}-fi  ■f*7"^4     7 rem     £$*#■  3£  !  . 

This  equation  is  good  along  the  line  9p=l80°  -  f  until  the  point 
©2=180°  -  7"  is  reached.  Bringing  Z  inside  the  parentheses  and 
recognizing  that:  £  s      gWjg  «?*+9') 

X  Of)         USr?    C  — — r       — 

6  *  o  fiP  * 


/.  *Vv)  r  ^ 


r  &Li  £ 


-  * /fe*i  r- 


rtOn 


It  may  be  noted  that  this  is  the  expression  for  magnitude  of  the 
first  overshoot  of  a  second  order  system  with  the  exception  of  the 
additive  constant   J^l1" » 

The  character  of  the  variation  of  the  magnitude  of  the  first 
peak  may  be  summarized  by  means  of  a  diagram  as  depicted  in  Figure lU. 
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The  coordinates  are  the  0z  and  0p  axes  used  for  the  fp   plots. 
The  line  0z=0p  has  been  drawn  in.  A  value  of  J*  =60°  (cosp"  =  .$) 
has  been  selected  arbitrarily  as  an  example.  The  value  of 
Mpt=l.l6  has  been  determined  from  the  second  order  equations  for 
placement  on  the  0z=0p  line.  From  Figure  6  values  of  Mpt  for  various 
0z  have  been  placed  along  the  0z  axis.  From  Figure  8(a)  values  of 
Mpt  for  various  0p  have  been  placed  along  the  0p  axis.  The  curve 
for  Mpt=1.0  has  been  copied  from  Figure  10a  (all  values  of  Mpt 
assuming  final  value  =1).  The  curves  outlining  the  regions  of  no 
overshoot  and  of  one  overshoot  only  have  also  been  drawn  in.  The 
lines  0z=l8O°-7"and  0p=l8O°-T  have  been  drawn  in  to  depict  the 
right  half  plane  (RHP)  and  left  half  plane  (LHP).  The  values  of 
overshoot  along  these  latter  lines  are  determined  by  the  equations 
derived  above. 

Approximate  curves  for  Mpt=l,l  and  Mpt=1.3  have  been  sketched  in 
to  show  the  form  of  the  curves.  This  form  is  that  the  curves  are 
convex  toward  the  0z=0p  line  and  all  curves  intersect  at  the  origin. 
The  value  of  Mpt  becomes  unbounded  for  values  close  to  the  0z=l8O- T 
line  and  approaches  zero  along  the  0p=l8O-T line.  Mpt  is  positive 
in  regions  I  and  IV  and  negative  in  regions  II  and  III.  It  is  a 
simple  matter  to  sketch  a  diagram  of  this  sort  for  any  value  of  T 
and  thereby  predict  the  action  of  the  system  both  in  time  and 
magnitude.  Similar  diagrams  for  the  first  undershoot  and  second 
overshoot  may  also  be  quickly  sketched. 
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No  curves  have  been  plotted  for  rise  time  (tc)  since  the 
computation  is  straightforward  and  the  information  already- 
plotted  gives  a  complete  picture  of  the  response  of  the 
system. 
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10.  A  New  Look  at  Some  Old  Concepts:  Dominance  and  Cancellation. 

It  has  been  stated  repeatedly  in  this  thesis  that  the  most 
common  method  of  treating  the  higher  order  systems  is  by  making 
a  second  order  approximation.  This  is  sometimes  done  by  forcing 
the  real  singularities  far  out  on  the  negative  real  axis  (dominance) 
or  by  cancelling  with  a  pole  and  a  zero  (cancellation).  Each  of 
these  methods  is  performed  in  a  rather  vague  manner j  i.e.,  how  far 
is  "far  out"  in  the  LHP  in  the  case  of  dominance  or  how  close  is 
"close  together"  in  the  case  of  cancellation.  Another  question  is: 
what  parameter  of  the  response  is  of  the  most  importance  in  making 
the  second  order  approximation?  Is  the  designer /analyst  interested 
in  time  of  response,  magnitude  of  fluctuations,  time  of  extremals, 
etc.?  It  is  possible,  by  making  use  of  the  methods  of  this  thesis, 
to  take  a  new  look  at  these  concepts  and  perhaps  make  them  a  little 
less  vague. 

Consider  the  case  of  a  third  order  system  with  a  troublesome 
pole  on  the  real  axis*  Using  the  concept  of  dominance  a  rule  of 
thumb  says  to  make  the  pole  ten  times  as  great  as  the  real  part  of 
the  conjugate  roots,  but  an  examination  of  Figures  7  through  9 
discloses  that  making  epsss30°  causes  a  variation  of  no  more  than 
10%   in  the  magnitudes  of  the  extremals  (the  magnitudes  are  less), 
no  more  than  1$%   in  the  times  of  the  extremals,  and  no  more  than 
30$  in  the  rise  time  (the  times  are  greater).  While  some  or  all 
of  these  variations  may  be  too  large,  two  points  are  demonstrated: 

(l)  If  the  time  parameters  of  the  response  are  the  important 
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factors  then  9p  must  be  made  smaller  than  if  the 
magnitude  parameters  are  the  important  factors. 
(2)  The  required  distance  of  the  pole  from  the  origin  varies 
with  the  value  of  Tf  .  For  example,  in  the  above  case, 
if  fl"»30°  then  it  would  only  be  necessary  to  make  the 
pole  twice  as  far  out  as  the  real  part  of  the  conjugate 
roots;  however  if  3*  =60°  the  pole  must  be  four  times  as 
far  out  as  the  real  part  of  the  conjugate  roots. 
It  is  fairly  certain  that  the  rule  of  thumb  value  of  "ten  times 
as  great"  will  make  the  pole  insignificant,  but  it  is  equally 
certain  that  it  is  not  necessary  to  be  so  drastic  in  all  cases. 
Utilizing  the  methods  described  here  and  armed  with  the  require- 
ments of  the  situation  at  hand,  the  designer  can  develop  the  most 
efficient  configuration. 

Consider  now  the  same  system  as  above,  but  assume  that  the 
troublesome  pole  is  to  be  removed  by  cancellation  with  a  zero. 
Examination  of  Figures  10  and  lU  show  immediately  that  the  require- 
ment of  how  close  the  zero  must  be  to  the  pole  to  effect 
cancellation  is  dependent  upon  the  location  of  the  pole  with 
respect  to  the  conjugate  roots,  value  of  the  angle  f,   and  the 
relative  importance  of  time  and  magnitude.  For  @p<70°  0z  can  be 
as  much  as  l£0  either  side  of  the  value  of  0p  and  the  time  of 
extremals  will  change  no  more  than  10$,  and  if  V  <   70°  the  change 
in  magnitudes  will  be  on  the  order  of  10$.  On  the  other  hand, 
if  6p>90°  the  value  of  6z  becomes  critical.  This  is  best  illustrated 
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by  examples 

Given:  T»60°;  9p=110°  (See  Figure  16) 

6z         Tp        Mpt 

11$°        2.8       1.958 

112°       3.0       1.353 

110°        3.1U      1.160  (exact  cancellation) 

108°       3.3       1.037 

105°       3.6       .908 

100°       U.55      .819 

100°      (no  extremal) 
For  a  value  of  P=.3  this  range  of  6z  represents  a  variation  of  Z 
from  .199  to  .U89.  It  is  readily  apparent  that  a  very  slight  error 
in  cancellation  can  cause  a  serious  deviation  from  second  order 
operation.  It  may  also  be  noted  an  error  of  6z>8p  causes  a 
smaller  change  in  Tp  and  a  greater  change  in  Mpt  than  an  error 
of  6z<8p#  Again,  it  becomes  easy  to  determine  the  margin  of  error 
allowed  to  effect  the  desired  response.  A  quick  glance  at  the 
charts  will  determine  for  the  designer  his  tolerance  and  which 
direction  to  weight  the  uncertainty. 
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11.  Conclusion. 

In  order  to  simplify  the  analysis  of  the  third  order  equation 
with  one  zero  the  time  response  equation  has  been  developed  in 
terms  of  the  angles  at  the  singularities:  ..  / j. 

This  equation  has  been  manipulated  to  provide  expressions  for  three 
time  response  parameters. 

(1)  Times  of  extremals  (  fan  )  ' 

Sin  (fifty) c'"/U"*/*=  <S*«  (utoJfi* +**-*/>) 

(2)  Magnitudes  of  extremals  ( X^/*) ) : 

(3)  Times  when  response  crosses  the  final  value  (  i£n   ): 

A  simple  graphical  method  has  been  developed  and  demonstrated  for 
solving  the  transcendental  equations  above. 

The  equations  above  have  been  utilized  to  analyze  four  types 
of  systems* 

(1)  Pure  second  order 

(2)  Second  order  with  one  zero 

(3)  Pure  third  order 
(U)  Third  order  with  one  zero 

Examples  have  been  worked  and  a  variety  of  curves  developed  to 
enable  the  designer  and/or  analyst  to  make  a  rapid  appraisal  of 
the  response  of  the  various  systems* 

An  analysis  has  been  presented  on  new  methods  for  achieving 
dominance  or  cancellation  and  evaluating  the  effects  of  errors  in 
such  attempts. 
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It  is  felt  that  a  new  insight  has  been  gained  into  the 
actions  of  the  third  order  system  with  one  zero.  With  the 
methods  and  charts  of  this  thesis  it  is  possible  to  predict 
immediately  the  form  of  the  response,  and  only  a  short  time  is 
required  to  compute  the  parameters  and  sketch  the  response  accurately, 
For  the  designer  interested  in  improving  the  response  of  a  second 
order  system  by  adding  real  singularities,  the  means  is  now 
available  for  doing  so  exactly  without  resorting  to  the  trial 
and  error  methods. 

There  is  a  possibility  that  the  methods  of  this  thesis  might 
be  extended  to  other  studies.  Some  that  have  occured  to  the 
author  are  listed  here  in  the  hope  that  others  may  be  persuaded 
to  attempt  to  extend  the  methods 

(1)  A  generalized  method  for  expressing  the  time  response 
of  higher  order  systems  in  terms  of  the  angles  at  the 
singularities. 

(2)  Combination  of  the  method  with  that  of  Mitrovic  to 
develop  generalized  curves  for  the  time  response  on 
the  Mitrovic  Plane. 

(3)  Extension  of  the  equations  of  this  thesis  to  develop 
charts  for  frequency  response  parameters. 

(U)  Combination  with  methods  such  as  those  of  Ross -Warren 

and  Pollack  to  develop  improvements  and/or  simplifications 
in  the  methods  of  compensation  calculations. 
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